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摘 要 在 多 排 叶 轮机 械 中 ， 动 静 叶 干涉 会 在 流 场 中 产生 强烈 的 非 定 常 效应 。 公开 文献 中 现 有 的 多 排 叶 轮机 械 气 动 
优化 设计 均 采 用 基于 混合 平面 法 的 定常 模型 ， 无 法 考虑 流动 的 非 定 常 特性 。 本 文采 用 课题 组 自行 开发 的 多 排 叶轮 机 械 
非 定 常 伴随 气动 优化 系统 ， 对 菜单 级 跨 音 压气 机 动 叶 进行 了 非 定常 气动 优化 设计 ， 优 化 过 程 中 采用 谐 波 平衡 法 高 效 求 
解 非 定 常 流 场 和 非 定 常 伴 随 场 .优化 目标 为 降低 压气 机 级 的 进出 口 箭 增 ， 同 时 限制 流量 和 压 比 不 低 于 原始 设计 ， 并 控 


制 叶片 最 大 厚度 的 变化 幅度 优化 结果 表明 ， 本 文 伴随 气动 优化 系统 能 够 高 效 地 对 多 排 叶轮 机 械 进行 非 定常 气动 优化 
设计 . 
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Unsteady Adjoint Aerodynamic Optimization of a Transonic 


Compressor Stage 


MA Can SU Xin-Rong YUAN Xin 
(Key Laboratory for Thermal Science and Power Engineering of Ministry of Education, Tsinghua University, 
Beijing 100084, China) 
Abstract In multi-row turbomachinery, the blade row interactions would impose strong unsteadiness 
in the flow field. However, to the authors’ best knowledge, all previous works on the aerodynamic 
optimization of multi-row turbomachinery are based on the steady mixing-plane model, which is not 
able to capture the unsteady flow features. In this work, the unsteady aerodynamic optimization 
of a transonic compressor stage is carried out with an in-house adjoint aerodynamic optimization 
system for multi-row turbomachinery. The harmonic balance method is used to efficiently solve the 
unsteady flow equations and the unsteady adjoint equations. The optimization target is to minimize 
the entropy generation from the stage inlet to the outlet, with no-decreasing mass flow rate and 


total pressure ratio. 


the maximum blade thickness. 


An additional geometric constraint is applied on the variation magnitude of 


Results show the efficiency and capability of the unsteady adjoint 


aerodynamic optimization System for multi-row turbomachinery. 
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基于 伴随 方法 的 梯度 类 优化 算法 在 叶轮 机 械 气 
动 优化 设计 中 得 到 了 广泛 应 用 号 3. 实际 的 叶轮 机 
械 优化 问题 中 设计 变量 一 般 较 多 ， 由 于 伴随 方法 的 
计算 基 与 设计 变 基 数目 呈 弱 相关 ， 采 用 伴随 方法 能 
够 高 效 地 求解 优化 所 需 的 梯度 信息 。 因 此 ， 现 有 的 
叶轮 机 械 梯度 类 优化 设计 中 大 部 分 采用 伴随 方法 . 
伴随 方法 分 为 连续 型 伴随 方法 和 离散 型 伴随 方法 两 
类 . 通常 ， 离 散 型 伴随 方法 的 梯度 计算 精度 高 于 连 
续 型 伴随 方法 四， 

在 多 排 叶轮 机 械 中 , 流动 本 质 上 是 非 定常 的 , 采 
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用 定常 模型 计算 会 存在 较 大 误差 尤其 是 当 叶 片 排 
数目 较 多 时 回 . 但 公开 文献 中 现 有 的 针对 多 排 叶 轮 
机 械 的 气动 优化 设计 均 采 用 定常 模型 计算 ， 无 法 考 
虑 流动 的 非 定常 特性 . 为 了 提高 优化 系统 的 可 靠 性 ， 
本 文采 用 课题 组 自行 开发 的 结构 化 有 限 体积 求解 
器 中 建立 了 基于 离散 型 伴随 方法 的 多 排 叶 轮机 械 
非 定常 气动 优化 设计 系统 ， 并 将 其 应 用 于 单 级 跨 音 
压气 机 二 维 叶 型 的 非 定 常 气动 优化 设计 . 为 了 提高 
求解 效率 ,， 非 定常 求解 采用 了 谐 波 平衡 法 门 , 优化 
结果 表明 ,本 文 伴随 气动 优化 系统 能 够 高 效 地 对 多 
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排 叶轮 机 械 进 行 非 定常 气动 优化 设计 . 
1 数值 方法 
1.1 流动 控制 方程 


考虑 有 限 体 积 方法 ,守恒 型 的 Navier-Stokes 方 
程 可 写成 以 下 形式 


其 中 UU 为 守恒 变量 
U = (p, pu, pv, pE)™ (2) 
尺 为 空间 残 差 项 , 由 对 流 项 下 。 和 生性 项 , 组 成 


1 
i Ee, (3) 
其 中 , Y 为 单元 体积 , n; 为 单元 面 法 向 , 5; 为 面积 ， 
采用 谐 波 平衡 法 时 ， 非 定常 流动 用 有 限 项 健 立 
叶 级 数 近似 ,求解 耗 时 相对 于 传统 时 间 推 进 法 能 降 
低 一 到 两 个 量 级 四. 设 非 定常 流动 角 频 率 为 w, 有 
oU 


Br 一 >》， ikw Ureit®t (4) 


k 
其 中 , i 为 虚数 单位 , e 为 自然 对 数 的 底数 . 在 谐 波 平 
衡 法 中 , 波 数 上 根据 流 场 中 非 定常 扰动 的 能 量 分 布 
选择 。 
在 一 个 周期 内 考虑 时 间 序 列 t;, j = 0,1,… ,2N 


记 
U*= [Uo, Di , UzN]! (5) 

则 
U*= EU, Ej= eK, (6) 
者 U*=E ikw E-'U* (7) 


ot 


应 用 谐 波 平衡 法 , 式 (1) 变 为 
DU*+R(U*)=0 (8) 


1.2 离散 型 伴随 方法 

离散 型 伴随 方法 从 离散 型 控制 方程 出 发 推导 相 
应 的 伴随 方程 . 假设 设计 变量 为 壁面 网 格 点 坐标 7， 
记 目 标 函 数 为 7(W, z), 其 中 W 为 原始 变量 


W = (p,u,v,p)T (9) 


目标 函数 对 设计 变量 的 梯度 可 写成 


dJ] 0J 0OJ OW 
Dm De 0 
对 于 定常 流动 , 流动 控制 方程 为 
R(W,7x)=0 (11) 
通过 引入 伴随 变量 4 满足 伴随 方程 
roR 07 
5W -Wm (2 
式 (10) 变 为 
dJ] 67J 10R 
i 0 


求解 伴随 方程 (12) 得 到 4 后 ,可 通过 式 (13) 高 效 
地 计算 目标 函数 对 设计 变量 的 梯度 . 

对 于 非 定常 流 动 ， 当 流动 方程 采用 谐 波 平衡 法 
离散 时 , 伴随 方程 在 式 (12) 的 基础 上 添加 物理 时 间 
项 , 形式 变 为 
oU” 


A*TD A*T 一 14 

Wt 于 
其 中 

4 入 一 [4o4 ,4ow] (15) 

W* = [Wo, Wi,:…, WoN]! (16) 


本 文中 ,流动 方程 与 伴随 方程 均 采 用 课题 组 独 
并 开发 的 结构 化 有 限 体 积 代 码 求解 . 对 于 流动 方程 ， 
对 流 项 计算 使 用 三 阶 WENO 迎风 格式 加, 务 性 项 计 
算 采 用 二 阶 中 心 格 式 , 潮流 模型 使 用 一 方程 Spalart- 
Allmaras 模型 . 伴随 方程 推导 过 程 中 未 考虑 清流 变 
量 对 庙 流 务 度 的 影响 ， 即 采用 定 涡 务 性 (Constant 
Eddy-Viscosity) 近似 ， 本 文 伴随 求解 器 对 黏 性 流动 
的 精度 已 经 过 验证 00 . 
1.3 优化 系统 

本 文 建立 的 优化 系统 如 图 1 所 示 , 在 每 一 个 优化 
循环 中 需要 求解 一 次 非 定 常 流动 方程 ,并 根据 目标 
函数 和 流 场 约束 条 件 的 个 数 求解 数 次 非 定 常 伴随 方 
程 。 本文 优化 算法 采用 商业 代码 SNOPTHH 中 提供 
的 SQP (Sequential Quadratic Programming) 算法 . 

在 本 文 优化 系统 中 , 参数 化 几何 建 模 采 用 FFD 
(Free-Form Deformation) 方法 站。 以 二 维 问 题 为 
例 ,，FFD 方法 定义 了 一 个 矩形 控制 体内 所 有 点 的 位 
置 变化 .该 方法 数学 形式 简单 ， 能 够 很 方便 地 得 到 
控制 点 坐标 与 几何 体 坐 标 之 间 的 梯度 信息 ， 特 别 适 
用 于 梯度 类 优化 方法 . 图 2 给 出 了 通过 FFD 方法 修 
改 压 气 机 尾 缘 几 何 形 状 的 示意 图 . 
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图 1 优化 系统 示意 图 


Fig. 1 Optimization System flow chart 


i 日 


图 2 优化 系统 示意 图 


Fig. 2 Optimization System flow chart 


几何 体形 状 发 生 改 变 后 ， 需 要 对 计算 域内 部 的 
网 格 点 坐标 进行 调整 . 本 文采 用 距离 反比 加 权 插 值 
方法 调整 计算 域内 部 的 网 格 点 坐标 ,任意 内 部 网 格 
点 的 位 移 5d 为 


0d = > Ga Da) (17) 
其 中 , 5d; 为 几何 体 表 面 及 网 格 边 界 点 的 位 移 , mr; 为 
内 部 网 格 点 到 前 述 网 格 点 的 距离 , p 为 插值 阶 数 . 网 
格调 整 过 程 中 进出 口 、 周 期 性 边界 保持 固定 ， 网 格 
块 之 间 的 连接 边界 随 内 部 网 格 点 移动 。 
2 优化 结果 

气动 优化 设计 对 和 象 为 菜单 级 跨 音 压气 机 [1 引 中 
截面 动 叶 的 型 线 , 其 几何 及 气动 参数 如 表 1 所 示 . 优 
化 采用 的 网 格 单元 数 为 1.0 万 , 原始 网 格 二 维 截 面 如 
图 3 所 示 . 优化 选取 的 工 况 点 为 最 高 效率 工 况 ， 原 
始 设计 压 比 为 1.28， 等 箭 效率 为 93.49%. 


表 1 压气 机 二 维 截面 参数 
Table 1 _ Parameters of the 2D slice of the 


compressor stage 


几何 及 气动 参数 动 叶 静 叶 
弦 长 c/mm 53.0 29.5 
栅 距 p/mm 32.8 15.3 
雷诺 数 Rec 4.9x105 3.4x105 

出 口 马赫 数 Ma 0.51 0.34 

平 动 速 度 v/(m/s) 229 0 
求解 基 频 f/Hz 1.5x104 7.0x103 
折合 频率 Foa 4.5 yi 


几何 造型 过 程 中 固定 前 缘 和 尾 缘 ， 只 对 中 间 部 
分 叶 型 进行 参数 化 造型 ， 从 而 保证 叶片 安装 角 、 折 


转角 和 弦 长 不 变 . FFD 控制 体 的 阶 数 为 8x8,， 如 图 
4 所 示 . 为 保证 叶 型 在 交界 点 处 一 阶 连续 , 靠近 控制 
体 边 缘 的 两 排 控制 点 在 优化 过 程 中 保持 固定 ， 实 际 
活动 的 控制 点 有 32 个 , 共 64 个 设计 变量 . 
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图 3 单 级 压气 机 网 格 (每 2 点 显示 1 点 ) 


Fig. 3 Mesh of the compressor stage (shown every two points) 


图 4 压气 机 动 叶 FFD 控制 体 
Fig. 4 FFD control grid of the compressor rotor 


非 定常 优化 的 目标 函数 为 计算 域 出 口 相对 于 进 
口 的 时 均 烂 增 


1 /™ 
As = 元 人 Sout UnsS; > punSidt— 
0 了 Sout 


>》 punssS; (18) 
| “om dt 
人 Jo >》 punS 

Sin 


其 中 了 为 周期 ， 下 标 s 和 + 分别 代表 静 叶 和 动 叶 . 
优化 目标 是 降低 时 均 炉 增 , 并 对 流量 、 压 比 进行 限 
制 , 要 求 优化 后 流量 、 压 比 均 不 降低 , 同时 限制 叶片 
最 大 厚度 变化 幅度 不 超过 0.1%， 
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非 定常 优化 过 程 中 参数 变化 历史 如 图 5 所 示 。 
优化 过 程 中 进行 了 179 次 流动 方程 求解 ，72 次 伴随 
方程 求解 ,求解 时 只 考虑 每 排 叶片 中 的 基 频 . 在 In- 
tel Xeon X5690 CPU 上 使 用 6 个 进程 并 行 , 耗 时 54 
h, 优化 后 炳 增 降 低 了 11.47%, 流量 、 压 比 均 未 降低 ， 
叶片 最 大 厚度 减 小 0.1%. 


19 20 


10 
优化 步 数 
图 5 非 定常 优化 参数 变化 历史 


Fig. 5 History of the unsteady optimization 


作为 对 比 ， 同 时 采用 混合 平面 法 进行 了 定常 优 
化 ,目标 函数 为 进出 口 炉 增 , 约束 条 件 与 非 定常 优 
化 相同 . 优化 耗 时 2.6 h,， 非 定常 优化 的 耗 时 约 为 定 
常 优化 的 21 倍 . 定常 优化 与 非 定常 优化 前 后 的 压气 
机 级 性 能 参数 变化 对 比如 表 2 所 示 . 两 组 优化 结果 
相差 不 大 , 等 炉 效率 仅 相 差 0.01%. 由 于 本 文 算 例 中 
静 叶 对 上 游 动 叶 的 影响 较 小 ， 动 叶 中 的 流 场 非 定常 
效应 很 弱 , 因此 两 组 优化 结果 相似 是 符合 预期 的 . 


表 2 优化 后 性 能 参数 变化 


Table 2 Performance variations after 


optimization 
定常 优化 非 定常 优化 
灶 增 11.54%] 11.47%} 
流量 0.62%1 0.71%1 
压 比 0.05%T 0.06%1 
等 粹 效率 。& 0.72%1 0.71%1 


优化 前 后 动 叶 叶 型 及 叶片 表面 时 均 压 力 系数 分 
布 如 图 6 所 示 . 压力 系数 定义 为 


Cp = (p — pin)/ (Puin — pin) 


非 定常 优化 与 定常 优化 结果 的 叶 型 和 表面 压力 
分 布 都 基本 一 致 , 优化 后 叶片 在 5% 到 40% 轴 向 粥 长 
之 间 负 荷 明 显 增加 ， 而 在 40% 到 80% 轴 向 弦 长 之 间 
负荷 明显 降低 . 同时 , 叶片 前 半 部 分 厚度 有 所 增 大 ， 
后 半 部 分 厚度 有 所 减 小 , 从 叶片 强度 的 角度 考虑 , 这 


一 厚度 变化 有 利 的 . 最 大 厚度 位 置 从 约 58% 轴 向 弦 
长 处 向 上 游 移动 到 约 48% 轴 向 弦 长 处 . 


原始 设计 
定常 优化 
a 非 定常 优化 


X/ Cax 


图 6 优化 前 后 动 叶 叶 型 及 表面 压力 分 布 
Fig. 6 Rotor blade profiles and surface pressure distributions 
before and after optimization 


非 定常 优化 前 后 的 瞬时 相对 马赫 数 和 炉 分 布 对 
比如 图 7 所 示 , 图 中 结果 采用 谐 波 平衡 法 计算 得 到 ， 
每 排 叶 片 中 计算 的 谐 波 数 为 5. 原始 设计 在 动 叶 吸力 
面 近 前 缘 附 近 存 在 局 部 超 音 区 域 ,峰值 马赫 数 达 到 
1.4, 同时 炉 值 较 大 , 最 高 达到 16 J kg KK- 优化 
后 该 区 域 峰值 马赫 数 降 低 到 1.1, 同时 最 高 粹 值 降低 
到 8Jkg 天- ,损失 明显 降低 . 另 一 方面 , 优化 后 动 
叶 尾 迹 也 明显 变 窗 , 同时 尾 迹 区 内 的 粹 值 明 显 降 低 . 


(a) 原始 设计 
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(b) 非 定 常 最 优 设计 
图 7 瞬时 相对 马赫 数 、 焙 分 布 


Fig. 7 Instantaneous relative Mach number and entropy 


distributions: (a) original design; (b) unsteady optimal 


design 


3 结 论 


本 文 以 课题 组 自行 开发 的 结构 化 有 限 体积 求解 
器 为 基础 ， 拱 建 了 基于 离散 型 伴随 方法 的 多 排 叶 轮 
机 械 非 定常 气动 优化 设计 系统 ， 对 单 级 压气 机 动 叶 
进行 了 二 维 型 线 优 化 设计 . 优化 后 压气 机 级 炉 增 降 
低 了 11.47%, 等 烂 效 率 提高 了 0.71%, 同时 压 比 、 流 
量 不 低 于 原始 设计 , 叶片 最 大 厚度 变化 不 超过 0.1%. 
优化 后 动 叶 前 半 部 分 负荷 提高 ， 同 时 厚度 增 大 ; 后 
半 部 分 负荷 降低 ， 同 时 厚度 减 小 。 优 化 后 动 叶 吸力 
面 近 前 缘 附 近 的 局 部 超 音 区 域 马赫 数 及 信 值 明显 降 
低 ， 同 时 动 叶 尾 迹 也 明显 变 窄 ， 尾 迹 区 内 炉 值 明显 
降低 . 作为 对 比 ， 同 时 采用 混合 平面 法 进行 了 定常 
优化 . 由 于 本 文 算 例 中 静 叶 对 上 游 动 叶 的 影响 较 小 ， 
动 叶 中 的 流 场 非 定常 效应 很 弱 ， 因 此 两 组 优化 结果 
相似 , 非 定常 优化 的 耗 时 仅 为 定常 优化 的 21 倍 . 优 
化 结果 表明 当前 伴随 优化 系统 能 够 高 效 地 针对 叶轮 


机 械 内 部 非 定常 流动 进行 气动 优化 设计 . 
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